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1 Motivations. 

The way moist-air surface heat flux should be com¬ 
puted in atmospheric science is still a subject of debate. 
It is explained in Montgomery (1948, M48), Businger 
(1982, B82), more recently Ambaum (2010, §3.7), that 
uncertainty exists concerning the proper formulation 
of surface heat fluxes, namely the sum of “sensible” 
and “latent” heat fluxes, and in fact concerning these 
two fluxes if they are considered as separate fluxes. 

It is shown in M48 that eddy flux of moist-air en¬ 
ergy must be defined as the eddy transfer of moist-air 
specific enthalpy h = Ci-at+pjp = eint + R F, where eint 
is the internal energy. However, the way the moist-air 
specific enthalpy h is computed in M48 still depends on 
some arbitrary assumptions concerning reference value 
of dry-air or water-vapour enthalpies, which are set in 
M48 to arbitrary conventional values at a finite refer¬ 
ence temperatures different from 0 K. 

Consequences of these arbitrary assumptions are 
studied at length in B82, though without succeeding in 
computing the moist-air enthalpy in an absolute way. 

Issues addressed in M48 and B82 can be overcome by 
using the specific thermal enthalpies derived in Mar¬ 
quet (2015, M15) for N 2 , O 2 and H 2 O, namely for the 
main components of moist air. 

In this article this approach is taken to show that 
Third-law based values of moist-air enthalpy fluxes is 
the sum of two terms. These two terms are similar 
to what is called “sensible” and “latent” heat fluxes 
in existing surface energy budget equation, but a new 
kind of “latent heat” is emerging in the definition of 
the moist-air enthalpy flux. Some impacts of this new 
“latent heat” flux are described in this brief version of 
a paper to be submitted to the QJRMS. 

2 The energy budget equation. 


the three local equations for ek, (j) and eint, yielding 
d 

^ [ P ( ^int + ^k)] = —V .[p{h-\-(l)-\-ek)U]-\-pq^ 

where g is a notation for impacts of radiation and other 
local sinks and sources of energy. The specific enthalpy 
h appears in the divergence term because the sum of 
— U.V(p) and —p V.U in equations for ek and eint, 
respectively, is equal to —V. (pU) = —V. [p{p/p) U], 
thus leading to the local definition eint + p/p = h. 



Figure 1: A column of atmosphere above Earth’s surfaee. 

The total energy budget equation is then computed 
by integrating this local equation over an infinite ver¬ 
tical column of atmosphere (see Figj^, leading to 

~~Qf~ ~ + 0 + + Q 

where lateral fluxes have been neglected. The term 
((...)) denotes the surface average oi p [h -\- (j) -\- e^) re, 
and the integral of p q ” is written as Q. 

The exchange of total energy between the column 
and the surface thus depends on { p {h -\- (/)) w ) which 
can be rewritten by using Reynolds average (...) and 
perturbation (...)' terms, leading to the local turbulent 
fluxes of specific enthalpy and potential energy 

Eh = {pwyE ~p w'h' , ( 1 ) 

= {p w)' (j)' ^ ~p w' (j)' . ( 2 ) 


Only three kinds of specific energies can be defined 
in atmosphere if nuclear or chemical reactions are not 
considered: i) the kinetic energy Ck = {u^ + v‘^-\-w‘^)/2] 
ii) the potential energy cj) — g z ] hi) the moist-air 
thermal internal energy eint = h — pjp or enthalpy 
h — Cint+p /p which are both associated with the First- 
Law of Thermodynamics. 

The total energy equation of a unit mass of moist air 
is computed for the sum etot = ^k + 0 + eint by adding 


Conclusions are the same as in M48 or B42: from 0, 
it is needed to know specific values of moist-air thermal 
enthalpy /i, in order to compute w' h '. 

3 The specific thermal enthalpy. 

The moist-air enthalpy is equal to the weighted average 
of individual (perfect gas) values for dry air, water 
vapour, liquid water and ice species, leading to /i = 
Qd hd -\-qvhv-\- qi hi -\-qihi. This sum can be computed 





















according to M15 (with different algebra), leading to 
h — /iref “ 1 “ ^pd ^ “ 1 “ Q.t -^vap Ql -^sub Qi 5 (3) 
where Lsuh{T) = hy{T) - hi{T) , 

Lh{T) = K{T) - hd{T), 

L..p{T) = K{T) - hi{T) , 
where qt = + qi is the total water content. 

The atent heats are in fact “differences in enthalpies”. 
They only depends on temperature and on some ref¬ 
erence values, with for instance Lh{T) = Lh{Tr) + 
(c-pv ^pd) ) 9.nd = {hv)r {h,d)r- 

Issues reported in M48 an B82 can thus be under¬ 
stood by a need to know {hd)r and {hy)r in order to 
compute Lh{Tr)^ then L^iT)^ and finally “/i” via (§. 
Dry-air and water-vapour reference thermal enthalpies 
are computed in M15 at 0 C: {h(i)r ~ 530 kJ kg“^ 
and {hy)r ^ 3133 kJ kg“^. The moist-air reference 
enthalpy is thus equal to = (^d)r ~ Cpd Ty ~ 
256 kJ kg“^. It is a true constant (whatever Ty may 
be in atmospheric range of T with constant Cp^) and 
it does not impact on gradient or flux computations. 



Figure 2: Comparison of Lh{T) withLyg,p{T) andLsuh{T). 
Unit are in kJ kg~^ for latent heats, in K for T. 


Changes of Tgub? and Tyap with absolute temper¬ 
ature are compared in FigJ^ The dashed straight line 
represents Lh{T). It is continuous at 0 C and is in 
between solid lines representing Lsub(T) and Lvap(T). 

4 The moist-air enthalpy flux. 

The moist-air thermal enthalpy fluxes ^ ~p w' N 
can be computed with h defined by (§, yielding 
Fh ^ Cp Ft + Lh Fy 

T — Z/vap) Fi — (Tgub “ ^/i) Fi , (4) 
where moist value of Cp is considered and where Fp ^ 
~p w' T', Fy ~p w' q'^ , Fi q[ and Fi g' 

are turbulent fluxes of temperature and water species. 

Clearly, — Tyap and Tgub ~ represent about 
4 to 8 % of Tvap? and about —8 to —10 % of Tgub- 
The second line of can thus be neglected with an 
accuracy better than 10 % and F^ ^ Cp Fp F Fy. 


However, if liquid water of ice contents exist, it is easy 
to avoid any approximation by taking into account the 
second line and the small terms depending on F) or F^. 

5 Numerical evaluations. 

From FigJ^and for Earth’s surface temperatures from 
0 to 30 C, Lh is about 6 % larger than Lyap in average. 
The change of Lvap(T) by Lh{T) should thus lead to 
an increase of about 6 % for existing surface latent heat 
fluxes which are about 100 W in average. Impact 
of using Lh{T) on global energy budget of atmosphere 
should thus be of the order of -h6 W in average. 
But impact larger than -h50 W m“^ could exist locally. 

These values have been confirmed by using short- 
range forecasts of ARPEGE NWP model, and in par¬ 
ticular for convective or frontal regions. 

6 Conclusions. 

It has been shown that computations of budget of to¬ 
tal energy imply computations of turbulent fluxes of 
specific enthalpy, and that it is possible to achieve the 
program started in M48 and B82 by computing these 
turbulent fluxes of enthalpy Q by using reference val¬ 
ues of enthalpies derived in Ml5. 

The surface enthalpy flux Q may replace what are 
commonly called “latent” and “sensible” heat fluxes, 
with the use of a new “latent heat” Lh{T) which is the 
difference in enthalpies of dry air and water vapour. 
Impacts of using Lh{T) is of the order of -h6 W 
in average, and more than -h50 W locally. 

Lh{T) represents physical processes like evapora¬ 
tions over oceans where water vapour enters atmo¬ 
spheric parcels of moist air, this meaning a decrease in 
qg at the expense of an increase of qy, in other words 
a replacement of dry air by water vapour. 

A striking feature imposed by the definition of the 
moist-air enthalpy is that Lh{T) is continuous at 
0 C. This is in direct contradiction to the usual def¬ 
initions of latent heat fluxes, assumed to be equal to 
^vap Fy over liquid water and Fgub Fy over ice. 
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